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Angular di-hadron correlations reveal novel structures in central Au-|-Au 
collisions at ^/Snn = 200 GeV. One of them, known as the ridge, is elon- 
gated in pseudo rapidity and peaks on the near side (A0 « 0). Investigating 
the origin of the ridge structure helps to understand the hot dense matter 
that is created in ultra relativistic heavy ion collisions. Results showing the 
(Pt) dependence of the ridge structure are presented. Evidence for possible 
jet and non-jet contributions to the ridge structure will be discussed. 

1. Introduction 

Recent triggered di-hadron correlation studies by STAR report a ridge 
structure in two dimensions (A?7,A(/))p]. Based on their studies, the near 
side was assumed to consist of two independent structures, a jet-like struc- 
ture and the ridge. A complementary approach, presented here, is carried 
out by using all possible charged particle pairs. The approach does not 
require a trigger particle but with appropriate kinematic cuts reproduces 
qualitatively similar correlation structure to that of triggered analysis. In 
this approach, two empirical 2d fit functional models are used to extract 
ridge properties. Possible model dependent quantitative evidence for jet 
and non-jet phenomena contribution to the near side structure can be ex- 
tracted from the fits. 



2. Data and analysis 

The data used in this analysis were collected during run 4, year 2004, 
from the STAR detector at Relativistic Heavy Ion Collider (RHIC) Brookhaven 
National Lab (BNL), Long Island, New York. 32M Au-|-Au collisions at 
\^Snn = 200 GeV were analyzed. Charged tracks reconstructed using the 
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Time Projection Chamber (TPC) with 2tt azimuthal coverage and \ri\ < 1 
in pseudo rapidity were selected. An earher centrahty evolution analysis [2] 
was based on particles in the full transverse momentum range above 0.15 
GeV/c. For the study presented here the lower threshold of the transverse 
momentum was raised for both particles. The selected event vertexes are 
within lb 25cm of the detector center and primary tracks are selected to 
come from within a distance of closet approach of 3cm to the event vertex. 

The correlation function for this analysis is a construct of Pearsons cor- 
relation coefficient. Mathematically it can be derived as[2j: 

Ap 

Psib Pmix 
\fPrej y/Pref 

Psib is the normalized charged particle pair density within a single event, 
which carries both correlated and uncorrelated pairs. The uncorrelated 
background is subtracted by mixed event pairs. The resulting uncorrelated 
density is denoted as Pmix • The difference is the number of correlated pairs 
which is normalized by the square root of mixed pair density to yield the 
correlated particle pairs per final state charged particle as the equivalent cor- 
relation measure to the Pearsons correlation coefficient. In order to avoid 
possible artifacts due to pseudo rapidity acceptance, the correlation mea- 
sure is carried out in sub bins of the z-vertex. The event mixing artifacts 
were eliminated by mixing events of a particular centrality within a multi- 
plicity window of 50 charged tracks. To correct for pair loss a pair cut is 
implemented on both mixed and same event pairs. The effects due to pileup 
events were removed by one of the standard pileup eliminating procedures 
in STAR[7J. Finally the conversion electron positron pair background is 
reduced by using a 1.5cr dE/dX cut on electrons in momentum ranges, 0.2 
<PT <0.45 GeV/c and 0.7 <pr <0.8 GeV/c. 

3. Fit functions 

The initial empirical fit function carries three mathematical model com- 
ponents that lead to seven free fit parameters. 

F = Co + cicos(A(/.) + C2exp(-O.5*((A0/c3)2 + {Ari/af) + C5exp(- 
O.5*((A0/C6)2 + (Ar?/C6)2)) 

The first component of the fit is an offset parameter followed by a 
cos(A(/>) term to extract the away side structure as seen in Fig. [TJ The as- 
sumption is that this structure corresponds to correlations due to a recoiling 
jet|3j. The two remaining terms are introduced based on {px) evolution of 
data (see Fig. [T]). The asymmetric 2d Gaussian attempts to address the long 
range correlation properties excluding the peak structure around (0,0) that 
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Fig. 1. Evolution of the raw normalized correlation structure for selected {px) bins 
at 0-1% centrality bin. The lower momentum threshold of both particles have been 
increased in generating the spectrum of correlation plots. 



appears as (jpt) increases. This structure is described via the symmetric 2d 
Gaussian component of the fit. 

Modifications to the model were introduced after indications of higher 
order harmonics were first predicted [H [5] and then observed (see Fig. [2]) . 

F = Co + cicos(A(/)) + C2Cos(2A(;/)) + C3Cos(3A(/)) + C4Cos(4A(/)) + 
C5Cos(5A(/.) + C6exp(-0.5*((A(/)/c7)2 + (/\r^/csf)) 

The resulting fit carries higher order harmonics up to 5th order and an 
asymmetric 2d Gaussian. It is important to note that in such a Fourier 
expansion, the first order term serves as a momentum conserving term 
(correlation due to recoil jet). In essence, the new model help in further 
constraining the knowledge gathered from the initial model study. 
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Fig. 2. Correlations at high pr very central data show evidence for the existence of 
higher order harmonics. 



4. Results 

The focus of comparing the initial model study to the study which in- 
cludes fit components for higher harmonics based on initial energy density 
fluctuations is to determine whether any information can be extracted from 
the asymmetric 2d Gaussian fit to the away side structure. 
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4.1- Initial model study 



Resulting asymmetric 2d Gaussian parameters are reported in Fig. [31 
The long correlation strength approaches zero at higher {pt) as indicated 
by the asymmetric 2d Gaussian amplitude and volume parameters. The 
Ar] and Acp width evolution clearly indicate the asymmetry shown in data. 
Also the At] width suggests that within the STAR acceptance the ridge is 
flat whereas the (j) width shows a very smooth narrowing. An important 
observation is that the (pt) of the ridge {{px) = 0.72 GeV/c) is greater 
than that of the inclusive spectrum {{pt) ^ 0.42 GeV/c [8j). The hardness 
of the ridge spectrum is an indication of possible Jet contributions to the 
ridge. 



I- Higher orders NOT included 
I- Higher orders Included 



: STAR preliminary 



<p^> [GeV/c] 



(a) Amplitude 



(b) Ar] width 



— AuAu 200 GeV 
^ CuCu 200 GeV 
PP 200 GeV 



- STAR preHmlnary 



0.5 1 1.5 2 2.5 3 3.5 4 



<p^> [GeV/c] 



(c) A(t) width 



STAR preliminary 
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Fig. 3. Asymmetric and Symmetric 2d Gaussian parameter evolution as a function 
of (pt)- 



On the other hand, the symmetric 2d Gaussian component reveals im- 
portant information about an unmodified jet emerging at higher momenta 
(Fig. [3]). The width is comparable to the jet width in PP collisions. In order 
to further constrain the above conclusion, amplitude and yield comparisons 
need to be made. 
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(a) ^ ratio (b) u„ scaling at — 10% 

Fig. 4. Comparison to predicted higher harmonic ratios[5] and scahng relations[3] 



4-2. Higher harmonics study 



Recent advancements in heavy ion collision models, predict significant 
contributions from higher order harmonics to the observed correlation spec- 
tra due to initial energy density fluctuations [H [5]. STAR data also sup- 
port the predicted observation in high pT, very central AuAu 200 GeV 
data(Fig. [2]). The ^ ratio predictions were tested and reasonable agree- 
ment to theory is evident (see Fig. [3^). However it is to be noted at low 
Pt, deviations from the trends are observed for Npart <150. One should 
note, thought, that the data is extracted fitting 2d di-hadron correlations, 
whereas, theory uses single particle spectra and event plane calculations. 

A predicted hydro scaling relation [6] using extracted higher order Fourier 
coefficients is studied, and reasonable agreement to theory has been observed 
(Fig. |3)3).In summary, given the different approaches used by theory pre- 
dictions tested in this proceedings and experiment, the data trends seem 
support the hydro based theory predictions to a reasonable extent. 

The remaining structure on the near side, after subtracting the Fourier 
components still shows a Ar]—A(f> asymmetry (see Fig. [3l red data markers). 
This {pt) evolution indicates that this Ar] elongation is visible up to 2.7 
GeV/c and is strongest at 0.9 GeV/c. The structure thus suggests possible 
jet modification at low (pt) which evolves to an unmodified jet at high pT- 
The amplitude and volume of the asymmetric structure scales lineary with 
Glauber scaling as a function of centrality, which hints at a jet origin. The 
rise of amplitude at high (px) is expected due to an increase in per charge 
particle pair yield in jets. Further studies will be carried out using the 
comparison to p+p collisions. 
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5. Summary and discussion 

Un-triggered 2d di-hadron correlation studies reproduce qualitatively 
similar results to that of a triggered analysis. The observed near side cor- 
relation is modeled via an empirical fit function which extracts short and 
long range structure properties. The {pt) evolution of the extracted pa- 
rameters suggests a possible correlation between jets and the observed long 
range correlation. Further constraint to data has been imposed modifying 
the empirical model to incorporate higher order Fourier model components 
cos(nA^), n = 1- 5, based on theoretical predictions and evidence in the 
data. The extracted higher order component strengths show reasonable 
agreement to predicted hydrodynamical trends. The remainder on the near 
side still reveals an asymmetric 2d Gaussian suggestive of possible modified 
jet phenomena. A Glauber linear scaling was successfully applied to the 
un-triggered centrality evolution (see Fig. [5|) which indicates a jet origin for 
the observed near side Gaussian. 
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Fig. 5. Comparison of near side asymmetric 2D Gaussian amplitude to Glauber 
linear scaling as a function of centrality. 
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